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[Text] Abstract 
A parallel matrix multiplier based on only a few imaging systems is presented. 


Many computation problems in science and engineering eventually become 
algebraic operations. Hence, the development of special devices to solve 
complicated algebraic operations has always been desirable. 


This paper presents a parallel matrix multiplier which involves only a few 
imaging systems. It is based on the corresponding decomposition of optical 
parallel vector-—vector out-product multiplier and matrix multiplier. It can be 
used to solve real analog multiplication of matrices, analog two-dimensional 
transforms, digital precision matrix multiplication, and digital precision 
two-dimensional transforms. It takes full advantage of the parallel nature of 
optics in all operations. The multiplication of three matrices or a two- 
dimensional transform takes place within one clock cycle of the system. 


I. Structure of the Parallel Matrix Multiplier 
1.1 Multiplying Two Matrices 


Matrix multiplication can be decomposed into a number of vector inner-products 
and outer-products (or tensor products). 


Figure 1 shows a novel optical parallel matrix multiplier based on vector- 
vector out-products. The system employs n point light sources on the focal 
plane of its spherical lens system. In the figure, SL and CL are spherical and 
cylindrical lenses, respectively. Their focal lengths obey the following 
equations: 








fa, *fa,* fa, *fa, =f, 








fox, = 3f/4 





In the figure, A and B are 
spatial light modulators 





(SLMs) used to input’ the g 
matrices. They are called ° 
input encoding plates. The SL, CL, Cla Cle Cls c 
transmittance of each element Figure 1. Optical Matrix Multiplier Based on 
is proportional to the value N Vector-Vector Outer Products 


of the corresponding matrix 

element. Therefore, we have to 

assume that every matrix element is a positive real number. As for different 
situations involving numbers that are not positive real numbers, there are a 
variety of methods available. Even a complex matrix can be coded to become a 
real matrix before input. In the system shown in Figure 1, addition and 
multiplication are done in the same manner as with an inner-product matrix 
multiplier. Figure 1 shows the principle of the optics. Obviously, it is an 
orthogonal imaging nonspherical optical system. 


If the light sources are ideal point sources and all lenses are ideal lenses 
(without taking aberration into account and cylindrical lenses only form 
images in one direction), then matrix code plates A and B are placed as shown 
in Figure 1. If the point source intensities are I)6(X-X,)6(Y-Y)), 1,6(X—-X,) 
5(Y-Y,), 12,6 (X-X,)6(Y-Y2), --+-:-> I n-1)6 (X-Xp_- ) 6 (Y-Yp-}), and the transmittances 
of encoding plate A and B are (Ajj) and (Bjj), then it is apparent that there 
is the following intensity distribution on plane P: 


n-1 
C13 * a) A y~Byy Ty 
0 


If Ip = I, s ceeeee - I,; = I, we have: 
n-1 


Cyy = Cty/GI = > Ay, Bay" (1) 


=0 


This is the product of matrix A and matrix B. Of course, transmittance is 
always less than 1. In reality, it is necessary to multiply by a proportion- 
ality factor. As for the rigorous derivation of equation (1), it can be 
obtained by using a transfer function technique. 


).2 Multiplication of Three Matrices 


If we need to multiply three matrices, then the system shown in Figure 2 may 
be used. the difference is the light source. In Figure 2, n x n point sources 
(assuming matrix B is n x n) are used. if the three matrices to be multiplied 
are (Ajj), (Bjj), and (Cjj;), then A and C are placed in the manner shown in 
Figure 2. The light intensity of each light source is proportional to the 
matrix Bj; (i,j = 0, l, crreee , n-l), and in their relative positions are also 
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shown in Figure 2. It is apparent that the light intensity distribution on 
plane P is: 





n-1 














where the symbol "~" repre- 
sents "proportional to." This 
is the product of the three 
matrices. Please note that the 
order of the three matrices in 





equation (2) is different from SLiCLi Cle Cle Cl 

that in Figure 2. This Figure 2. Optical Three-Matrix Multiplier 
arrangement provides some Based on Vector-Vector Outer Product 
convenience for the 


multiplication of more matrices. 


II. All-Parallel Optical Multiplier Technique and Its Optical Two-Dimensional 
Transform 


Because light intensity is always positive, all matrix elements described in 
the previous section are positive. In the multiplication of matrices, some 
elements may be real or complex numbers. This requires us to take certain 
measures. 


This paper presents a technique which involves a combination of polarization 
encoding and light intensity encoding to realize optical two-dimensional 
transform. 


2.1 Expression and Operation of Intensity/Polarization Hybrid Encoding for 
Real Numbers 


(1) Expressing a real number 


To express a real number by intensity and polarization, the magnitude of light 
intensity can be used to express the absolute value and two different 
directions of polarization can be used to represent different signs. Based on 
this concept, an arbitrary real number may be expressed as follows: 


M = |M|[1 + sign (M)]/2(+1) + |M|[1 - sign(M)]/2(-1) 


where 
+1 MO 
sign(M) = 40 M=0 (3) 
-1 M<O 


Optically, e, and ey, are two independent orthogonal polarization directions. 
It is further assumed that e, corresponds to (+l) and e, to (-1). Then, the 
light beam representing the number M has the following two polarization 
components ({I,, Ty}: 





I, = |M|(1+sign(M)]/2 (4) 
I, = |M|{1-sign(M)]/2 
Thus , Mom Iyex + lyey 
and 
I,ly = 0 (5) 


Obviously, this expression is unique. In order to meet the need in practice, 
the two independent polarization components of the light beam expressing M can 
be expressed in a more general manner as: 


I, = |M|[1+sign(M)]/2+C (6) 
Iy = |Mj(1-sign(M)]/2+C 
C20 


In this case, the expression is not unique. 
(2) Arithmetic operation 


The multiplication of two real numbers can be expressed by the method shown in 
equations (4)and (5), i.e. 


M = Iyxex + Iyyey 
N = Iyxex + Iyyey 
Ty = |Mi(1+sign(M)]/2 
Iy, = |M|(1-sign(M) ] /2 
Iy, = |N|(1+sign(N)]/2 
Iy, = |N|(1-sign(N)]/2 
MN = |M||N|[1 + sign(M) sign(N)]/2e, 
+|M||N|[1 -— sign(M) sign(N) ]/2ey 
= |MN|[1 + sign(MN)]/2e, + |MN|[1 - sign(MN) ]/2e, (7) 


The multiplication of two real numbers optically expressed by equations (4) 
and (5) can be realized by shining two polarized light beams on two 
transmission plates. If the transmittances are ,, m2, then 


I = myn2lToxex + Ioyey | (Toxloy = 0) 


In order to express the sign of operation of the two numbers, it is possible 
to stick a layer of polarizer on the transmission plate. Transmittance and 
polarization angle together can express a number M with a sign. When M < 0, 
then a layer of polarizer can be used to rotate the direction of polarization 
of the incident light by 90°. When M20, no polarizer is required and the 
direction of polarization remains unchanged. In the meantime, n, = |M|. The 
intensity of the incident beam represents the value of N and its direction of 
polarization represents the sign of N. Then, the output is (M-N), which in 
equation (7) is the product of two real numbers. 

















If addition or subtraction is required, then it is necessary to use the method 
shown in equation (6). In this case, the expression is no longer unique. For 
example, the subtraction of two greater than 0 real numbers A and B, i.e., 

C = A — B, can be expressed as: 


C= Tax€x + Igyey 


In this case, subtraction becomes adding a positive value to the polarization 
component which represents the negative value. 


Obviously, the number must be transformed after addition or subtraction before 
multiplication can be carried out again. Nevertheless, it is very difficult to 
do so directly by optical means. Hence, it is usually done electrically. 


2.2 Two-Dimensional Transform of Intensity/Polarization Hybrid Code 


Any two-dimensional linear transform in the complex domain can be reduced to 
four independent real number transforms. Therefore, this section only deals 
with two-dimensional linear transforms in the real number domain. (This 
remains true unless specifically expressed otherwise. ) 


Any two-dimensional transform can be expressed as follows: 
E = UFV 


where F is the matrix to be transformed, U and V are the transformation core, 
and E is the end result. If Fjj is expressed by the method described earlier. 
|U;;| and |Vjj| are expressed in terms of encoding plate transmittance and sign 
(Uj;) and sign (V;j;) are expressed by the polarizer sticker method. The optical 
system to realize this transformation is shown in Figure 2. Its light source 
is: 


Fi; = |Fij|6(x-x;) 6 (y-yj) ([1 + sign(Fjj)]/2e, + [1 - sign(Fjj)]/2ey) 
The encoding plates are: 


Uij = [Uji | 6 (x-xi) 6 (y-yi) ([1 + sign(Ujj)]/2e, + [1 - sign(Uj;)]/2e,) 
Vij = Vig) 6 (xxi) 6 (y-yi) ([1 + sign(Vjj)]/2e, + [1 - sign(Vjj)]/2e,) 


For a specific type of two-dimensional transform, Uj; and Vjj remain unchanged. 
Hence, they can be constructed as fixed encoding plates. Thus, the system does 
not require an SLM as the input; it only needs two light arrays that can be 
modulated. If Fj; is always greater than 0, it is even simpler. It only needs 
one light source array, which is extremely convenient. 


In Figure 1, all we have to do is to subtract the two polarization components 
received correspondingly. This only requires equation (6), and the result is 
apparent. The subtraction of two polarization components can be completed on 
an electronic circuit. There is no need to do it optically. 








In conclusion, the principle of an optical parallel vector-vector outer- 
product multiplier is presented. It is very simple and can easily be expanded 
to larger dimensions. By using common optical components, it is possible to 
obtain the outer—product of 10% x 10‘ vectors (in analog). The system can be 
used to realize multiplication involving two and three (positive real number) 
matrices (in analog). Experimental results to verify this theory will be 
published separately. (See next issue.) 


On the basis of hardware research on an analog matrix multiplier, the 
algorithm is discussed in this paper. It expands the range of applications to 
the multiplication of real number matrices. Furthermore, a hybrid spatial 
intensity/polarization encoding method is used to realize the multiplication 
of two to three digital (mixed-base) matrices to improve the accuracy of 
operation. Compared to the two existing all-parallel matrix multipliers (based 
on inner-product and Fourier transform), this system has a distinct advantage 
in performing two-dimensional linear transform in the real number domain. It 
only needs a light array to represent the pattern. The transformation core is 
a number of fixed encoding plates. The result is obtained in a single step. It 
does not require an SLM. (This is because matrix B, in the multiplication of 
matrices A, B, and C, is expressed by a light array in this approach. In other 
approaches, matrix A is expressed by a light source instead.) Therefore, use 
of this triple matrix multiplier system should be promoted for two-dimensional 
transforms. 
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[Text] Abstract 


The experimental method and results of a parallel optical matrix multiplier 
are described. The system takes full advantage of optical parallelism and is 
capable of handling the multiplication of three matrices of a two-dimensional 
transform within one clock cycle. The sources of error and ways to eliminate 
them are discussed. Results of multiplications involving two and three 
matrices (analog operation) are given. The mean error of two-matrix multipli- 
cation is 1.2 percent and that of three-matrix multiplication is 1.23 percent. 


Currently, matrix operations are done serially in computers. Hence, the 
inherent parallelism of matrices is totally lost.!* Optical computing is 
extremely parallel in nature. Therefore, researchers are actively seeking 
optical means to realize parallel matrix operation. To date, there are many 
optical matrix multipliers which have been briefly discussed in reference [3]. 
A parallel matrix multiplier involving only a few imaging elements was 
introduced in reference [3]. Multiplication of real number matrices is 
realized based on optical parallel vector-vector outer product and the 
corresponding decomposition of matrix multiplication. This paper describes the 
experimental process and results of this parallel matrix multiplier. It makes 
full use of optical parallelism. The multiplication of three matrices or a 
two-dimensional transform takes only one clock cycle. 


Principle of the All-Parallel Matrix Multiplier 


Figure 1 shows the vector-vector outer product based optical parallel matrix 
multiplier presented in reference [3]. The system employs n x n point sources 
placed on the focal plane of a spherical lens system for illumination. 








In addition, it represents the matrix B. In the figure, SL and CL are 
spherical and cylindrical lenses, respectively. Their focal lengths follow the 
equations below: 


fot, = fa, = fo *fa=fifas* 3f/4 


The SLMs used to 
input the other 
two matrices are 
A and C. In our 
preliminary 
experiment, they 
were replaced by 
input coded 
plates. (In a 
two-dimensional 
transform, these 
two coded plates 
represent the 
core of trans- 
formation.) The 
optical trans- 
mittance of each 
element is pro- 
portional to the 
value of the 
corresponding i) 0) 
matrix element. 
To this end, we 
have to assume 
that all matrix 
elements are 
positive real numbers. When they are not positive real numbers, there are a 
number of schemes available. Even when they are complex numbers, a matrix can 
be coded to become a positive real number matrix‘ prior to processing. In the 
system shown in Figure 1, optical addition and multiplication are accomplished 
in the same manner as calculating the inner product of matrices. Figure 1(a) 
illustrates the principle of the optical system and Figure 1(b) shows the 
paths of light rays. Apparently, it is a nonspherical orthogonal imaging 
optical system. 


























Figure 1. (a) Optical Scheme of the Three-Matrix 
Multiplication System 
(b) Paths of the Light Rays 


If the light sources are ideal point sources and all lenses are ideal imaging 
lenses (without considering any aberrations, and the cylindrical lens forms an 
image only in one direction), then coded plates A and C are placed according 
to the way shown in Figure l(a). If the light intensities of the point sources 
are By,, the transmittance values for coded plates A and C are (Ajj) and (C;j), 
then the receiver plane should have the following intensity distribution: 


Diy = GEA, )BypCyyDyy = Djy/G = BAy)By yyy" (1) 














This is the multiplication of matrices A, B, and C. Of course, transmittance 
is always less than 1. In reality, it is always necessary to multiply by a 
proportionality constant. 


For the multiplication of two matrices, all we have to do is to change matrix 
B to a unity matrix, i.e., to change the light source to n point sources. 


Analysis of the Experimental All-Parallel Optical Matrix Multiplier System 


The major sources or error for the matrix multiplier described above include: 
(1) error due to diffraction of small elements on the coded plates, (2) error 
due to nonuniform optical angular distribution of different point sources, 
(3) error from background light, (4) error caused by distortion of coded 
plates, and (5) error due to deviation of the actual optical system from an 
ideal system. When the size of the elements of coded plates are large, (2), 
(3), and (4) have larger impact on the error. It is easy to reduce the last 
one to meet the accuracy requirement. 


1. Relation between size of optical elements and coded plates 


Let us assume that the point light sources are s, and s, from center to center 
and the elements on plates A and B are dj, dy, dx and d,, from center to 
center, respectively, and each element of receiver is a small block h, x hy in 
size. If the system is an ideal optical system, then based on optical wave 
theory, the resolution light is 


hy = 1.22A/NA, 
hy; = 1.22A/NA, 
Let hy; - hy; _ h, 
Then h, = hy = h>2h, 


It is obvious that 
Sx = Sy = dix wdax = diy = day = h 


The aperture of the system should be determined by the outer frame of A and B. 
Hence, there are specific requirements on minimum aperture for every lens. 


In order to avoid any overlap of images between neighboring elements on the 
receiver plane, it is necessary to have an opaque area between neighboring 
elements. The width of this area is determined by the resolution of each 
component of the system. Similarly, when each light source is not an ideal 
point source, the size of the light sources and the gap between them are 
determined by the resolution of the system and the element size on the 
encoding plates. Because optical wave theory can only provide a rough 
approximation with regard to the data (and it often is far off from reality), 
it will not be discussed in detail. Experimentally, the spacing between two 
neighboring elements was chosen to be R times the width of the element itself 
(when approaching the resolution limit of the system, R > 1; when far away 
from the resolution limit, R << l). 











2. Error due to nonuniform angular light distribution and ways to eliminate it 


The error due to nonuniform angular distribution of a light source is obvious. 
With Gaussian distribution, it can be expressed as: 


(6) = I,e-” 


When the light from a point source passes through the matrix multiplier and 
reaches A, it can be expressed as follows: 


I,y~e -ay*/f? 


y~-ftg 6= fé 


After penetrating plate A, it is 
I,y~e*” /*"f , (y,) 


Upon reaching B, it becomes 
Iyy~I,yo*"/® 


After penetrating B, it is 7, ~r,,e-«/ff,(x,) . On the receiver plane, we have 


I yiy1 * Ge LY N/LF (ys) fg (X,) 
where G is dependent upon the diaphragm of the encoding plate. 


In order to allow this distribution to meet the accuracy requirement 6, we 
have 


f/D = [18a/1n(6+1) }?/? 


where a is a parameter which defines the Gaussian distribution, § is the 
accuracy requirement, and D is the maximum width of the encoding plate, i.e., 
D = max (Dx, Dy). In practice, different light sources may be chosen based on 
the data. Of course, it is only feasible based on availability. 


3. Methods to eliminate error 


Spatial distribution multiplication error can be eliminated by a simple 
method. It involves adding encoding plates with Ajj = Cjj = 1 and Sjj = 1 onto 


matrix plates A and C. Use the receiver to measure pj, and determine 
min( pi,)/ Di; = Djj. Place an encoding plate with transmittance Dj; immediately 
in front of the receiver. This is a way to eliminate any multiplication error. 


Similarly, additive error may be eliminated electronically at the receiver. 
Using this method, error can be minimized to a very low level (below 1 per- 
cent), because most errors in the optical system can be eliminated by these 
two techniques. 
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System Adjustment, Experimental Procedure, and Results 


Using the system shown in Figure 1, experiments were done to multiply two 
matrices AioxsBsx10 and three matrices Uroxs, Vsx5> and Wsx10- 


1. System adjustment 


First, align the center of each lens and encoding plates. Adjust each lens so 
that it is perpendicular to the center line. Then, adjust the two sets of 
cylindrical lenses to make them orthogonal to each other. Based on the 
dimension and accuracy requirement of each optical element and encoding plate, 
the following parameters can be obtained: 


center deviation < (d-D)/4; 
vertical deviation < h/D/2 


where d is the actual aperture of the optical element, D = max(D,, Dy), and h 
is the minimum width of an encoding plate element. 


2. Experimental procedure 


In the experiment, f = 10 mm, each encoding plate element is approximately 
0.6 x 0.6 mm’, and the gap is 0.4 mm or 0.3 mm wide. Every cylindrical lens is 
mounted on a two-dimensionally adjustable rack. Each encoding plate is clamped 
onto a three—dimensionally adjustable rack. Experimentally, five point sources 
in Bj; were first used for illumination. After using encoding plates Ajj = 1 and 
Cjj = 1, adjustments were made to obtain uniform results on the receiver 
surface. Then, the real encoding plates were inserted to conduct the tests. 
The error of the display is 1.0 percent. 


3. Experimental results 


(1) Multiplication of two matrices: 
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a 


bs 

q 

" 
FN WD PF FP DY WD PE 
NWNPRP RP ND WY DP Pb 
WDNR rR DW DFP FP DY 
NrePrPN WN FP FD W 
KFrPN WN FP FD WY DN 
RPN WD PP DY WD PE 
NWN PRP PRP DW DP 
WNnNrP FN WD BS FD 





rPrRPNWNR RP DD WN 





NrPrP ND WDNR FD W 


(2) Multiplication of three matrices: 
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1134311343, 








Figure 2 [photo not reproduced] shows the result of multiplying two matrices. 
The picture was enlarged (and a correction plate as the one described earlier 
was inserted in the experiment). For convenience, each element is a plate with 
either 1.0 or 0.0 transmittance. Table 1 shows the optoelectronic measurement 
results for the situation shown in Figure 2 [photo not reproduced]. 
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Table 1. Results of Two-Matrix Multiplication (Detected by optoelectrical 





detector) 
8.01 2.02 0.98 1.01 1.98 3.08 2.04 0.99 1.03 2.02 
1.97 2.98 1.97 1.02 0.99 1.99 2.90 2.01 1.02 0.98 
0.99 2.02 3.03 2.01 1.01 0.98 2.08 8.04 2.08 0.99 
1.02 1.01 2.04 8.05 1.97 0.99 1.01 2.08 3.01 1.97 
1.98 0.98 1.01 2.03 2.97 1.99 0.99 1.01 1.99 3.02 
8.06 2.01 0.99 0.99 2.01 2.97 1.96 0.98 1.01 1.98 
2.02 3.01 2.03 1.02 1.01 1.99 2.96 2.09 0.98 0.99 
0.99 1.97 2.95 1.97 0.98 0.99 2.08 3.05 2.03 1.01 
1.00 0.99 2.01 2.95 2.04 1.01 1.02 2.08 3.03 2.04 


1.99 0.97 1.01 2.01 3.02 2.03 1.01 1.00 2.01 3.04 





Errorges"2% Error... 1.2% 
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[Text] Abstract 


Optical multistability and optical ternary logic gates have been realized 
using a composite ZnS interference bistable device. The device is an ideal 
element for studying the feasibility of optical computing because of its 
favorable characteristics such as ease of fabrication, low threshold power, 
and room-temperature operation. 


Interference-filter—based optical bistable devices are simple to fabricate and 
easy to operate. They may be used as logic elements in areas such as optical 
communications, optical computing, and optical information processing. 
Interference filter—based optical bistable devices have been successfully used 
in optical full adders, mode recognition,? digital conversion,’ and high-speed 
switching.‘ This paper reports the use of a composite interference filter 
bistable device to realize optical multistability and optical symmetric 
ternary logic gates. Compared to binary logic systems, a multi-valued logic 
system has higher operating speed and information storage density. Recently, 
Liu Shutian and M. Watanabe investigated optical multi-valued logic devices 
such as the Post ternary logic gate and symmetric ternary logic gate using 
optoelectronic mixers,*>® respectively. Their work has’ consice~able 
significance in improving the capability of optical computing. 


I. Optical Multistability 
A dielectric interference filter contains alternate coatings of high- 


refractive index material (such as ZnS or ZnSe) and low-refractive index 
material (such as Na,AlF, or MgF,). Each layer is one-quarter wavelength thick. 
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The film structure is (HL)™(nH)(LH)™. (HL)™" and (LH)™ are considered as two 
highly reflective mirrors on both sides and the entire filter is equivalent to 
a nonlinear F-P standard device. The ZnS interference filter is intrinsically 
a thermally activated optical bistable device. The bistability is caused by a 
shift of the center wavelength of the filter to the incident wavelength due to 
heat. The threshold power to produce bistability is dependent upon the 
incident angle. The larger the angle of incidence is, the higher the threshold 


power becomes. 
Ls 


Bs Li L; 








In this experiment, two interference 

filters were used in series to obtain Ce 
optical tristability. The experimental 
setup is shown in Figure 1. The output of 
the first interference filter is used as 
the input for the second filter. The fil- Figure 1. Experimental Setup 


ters have a center wavelength of 513.1 nm Ar*: argon laser; M: modulator; 
and 513.8 nm, a FWHM of 2.1 nm and 1.7 nn, BS: beam splitter; 

and minimum threshold power of 15 mW and L,~L,;: lenses; IF,, IF,: 

11 mW, respectively. Figure 2 [photo not interference filters; 
reproduced] shows a tristability curve D,, Dz: detectors 


obtained by varying the input light power. 





The tristability produced is in reality ‘ 
the superposition of two bistable 
effects. As the incident light intensity 
increases, the output light intensity ? 
also continues to rise. When the input 
light intensity, P;, is equal to that of 
the threshold, IF,, IF, is opened. At 
this time, although IF, is still at its 
lower state, the output light intensity, 
Py), still jumps considerably because of 
its relatively high background transmit- 
tance. As the input light intensity in- 
creases, when the light intensity pass- FP 
ing IF, reaches the threshold of IF2, IF; 

is opened and the output light intensity Figure 3. Truth Tables and Input- 
jumps up once again. This is how tri- Output Relation Curves for 
stability is generated. The shape of the "OR," "AND" Gates 

curve may be varied by changing the 

incident angles of IF,, and IF3. 
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II. Optical Multi-Valued Logic 


The most basic operations in symmetric ternary logic are "AND" and "OR." 
Figure 3 shows the truth tables and the input-output relation curves for these 
two logic gates. It indicates that OR and AND gates can be realized by the 
same optical path. It only requires changing the thresholds of the two 
bistable devices. Under conditions to meet the logic relation for the AND 
gate, it is possible to meet the logic relation for the OR gate by raising the 
threshold of the first bistable device and lowering that of the second device 


15 





at the same time. Experimental results of the AND and OR gates obtained using 
the optical system illustrated in Figure 1 are shown in Figure 4 [photo not 
reproduced]. Pulses of different amplitudes correspond to different input 
light intensities. P,, P,, and P2 are the discriminating intensities of 0, l, 
and 2, respectively. 

IF; 
The NOT gate [or "inverter"] is also a common 
logic gate. It may also be realized with two 
filters, as shown in Figure 5. The light 
reflected by the second filter is the output. —. 
Due to polarization P,, when the input is 0, l, IF, 
and 2, the output is 2, 1, and 0. This meets 
the logic relation for the NOT gate. P, 





ov 





III. Discussion 2 


Interference filters are used to study the 

feasibility of optical computing. Later, they 0 

can be replaced by high-speed switches. Because 

optical systems have a high degree of parallel- tS P, 
ism, high-processing-speed optical computing 
systems of practical value can be achieved even 
with medium-speed logic elements (1 ps~l1 ns). 
The use of medium—-speed logic elements is also 
compatible with available light sources and signal detectors. As progress in 
interference filter bistable devices advances, ZnSe filters operating in the 
infrared have been successfully developed. Such devices can be used with laser 
diodes to substantially reduce the volume of the system. 





Figure 5. Principle of NOT 
Operation 
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[Text] Abstract 


The optical parallel cache memory array and random access array are critical 
parts of a digital optical processor. A simple optical parallel cache memory 
array is presented. Furthermore, random access to the optical array has been 
achieved experimentally. 


I. Introduction 


Optical mathematical computation is a highly parallel spatial optical logic 
operation.! It allows the simultaneous transfer and processing of a vast amount 
of optical information through a high-density parallel optical network and 
two-dimensional optical switching array.* From the algorithm and system 
architecture of digital optical computing, the two-dimensional optical logic 
array is essential to such a parallel digital optical processor. The optical 
parallel cache memory array and parallel random access memory array are also 
vital to such a system.’ Because digital optical computing is highly parallel, 
in order to suit high-speed optical processing, the data, command, and program 
required for processing are first stored in the cache memory array for ready 
access. In optical information processing, there is also a need to store 
intermediate results in random access memory to be taken out in the next step.‘ 


There are numerous studies on optical bistable logic devices.°*:® In this work, 
a random access memory array consisting of two-dimensional optical bistable 
devices is presented. Details regarding writing, reading, and zero setting are 
described. Parallel optical random access of information was experimentally 
achieved. 
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II. Optical Bistable Memory 


In digital optical computing, an optical bistable device can be used as an 
optical logic switch, as well as optical cache memory array.’ As an optical 
logic device, its bistability curve should be adjusted to approach 0, as shown 
in Figure 1. As an optical memory device, it requires a more stable 
bistability range in order to implement basic functions such as _ random 
writing, storage, reading, and clearing. 


In digital optical 
Pout computing, the binary 
code of spatial logic 
° consists of a series bY 
of black and white c 
b patterns.!:* A bright 
—<— 

. Pin spot has a certain 
level of power. A 
dark spot has very a Pin 
low power, approach—- Figure 2. Output/Input 
ing 0. For the sake Characteristic 
of discussion, let us for the Memory 
assume a bright spot is "1" and a dark spot is 
"0." When a bistable device is used as a memory 
device, as shown in Figure 2, the power of bias light H, is first adjusted to 
point b in the bistable zone. When the input signal is "0," point b belongs to 
the lower state and the output signal is "0." When the input signal is "1," it 
jumps from point b to point c which is at a higher state and the output is 
"1." After the signal is removed, the output still remains at point b’. There- 
fore, the original optical signal "1" is stored in memory. If this information 
needs to be randomly accessible, an additional optical bistable switch is 
required. Its characteristic curve is shown in Figure 1. The bias light for 
this device is set at point b before the bistable zone. If one needs to take 
out the information at a certain point of the memory, the bias light must be 
applied to the corresponding point on the switch to read it out. 


Pout | 

















Figure 1. Output/Input 
Characteristic for 
Logical Switching 










































































” If the next piece of information needs to 
be stored at the same location, then the 
aisle data stored there must be cleared. This can 
bby 6 array ST be done by cutting off the bias light H, on 
== _ ¢ or the memory and the corresponding bias light 
Bs ES 1/12 H, on the switch to set 0. 
AAN ] 
RSS g3 cat BIE 
sie oan2 III. Optical Parallel Cache Memory Array 
NU 4 ty & array of aicroles Experiment 
SY 
Re orien Figure 3 shows the experimental setup for 


the optical system. The optical bistable 
device is the normal transmission-type 
optical bistable ZnS device®:® developed in 
our laboratory. Its bistability character- 
istic curve is shown in Figure 4. 


Figure 3. The Experimental Setup 
of 4 x 2 Array for Parallel 
Optical Cache Memory Array 
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The x and y axes are the input and output power, respec-— 
tively. In the optical path, the planar thin-film beam 
splitter splits the beam into a 4 x 4 array. The SSR 
composite prism and the 4 x 4 microlens array were all 
developed in our laboratory. 


The SSR composite lens can divide the 4 x 4 beam array 
into two arrays of identical intensity. Moreover, the 
two beam arrays can be displaced relative to each other 
and make beams from different points in space coincide 
with one another in a rigorous manner. In this experi- Figure 4. Bistable 
ment, the first and third column and the second and ¢yaracteristic of 
fourth column were made to overlap by the SSR prism to ywo;na1] Transmission 
form a 4 x 2 array, as shown in Figure 3. Using the half 27,6 IF pevice 

wave plate HWP,, the power of each beam in column P is 

adjusted to 9 mW. P, and P- are used as the bias light 

H, for the memory. In the S light array, the S,S, array is the signal light 
source. Each beam has a power of 3 mW. Arrays S,S, and P,P, are overlapped and 
passed through the microlens to enter the normal transmission optical bistable 
device OBD, at right incident angle. Array P3P, is used as the bias light H, for 
the second optical bistable device OBD,. Its optical power may be adjusted by 
HWP, . 





In the experiment, a 
moving screen M, was used 
to block off light beams 
in columns 3 and 4. 
Therefore, P,P, is the 
only bias light shining 
on OBD,. If the bias 
light is modulated sinu- 
soidally from 0 to 9 mW, Figure 5. The Process of Writing to, Reading From, 
the bistable curve of and Clearing the Memory in Optical Cache 
the 4 x 2 array output Memory Array 

of OBD, is as shown in 

Figure 5(a). From the photograph, when the bias light power is at 9 mW for 
every beam, every point of the 4 x 2 array is at its low state. When the input 
light signal is "0," the OBD, output signal is "0." When M, is used to block 
off the light array P,P, and M, is removed to let the light array S,S, with 
optical signal "1" through, every point of the 4 x 2 light input array for OBD, 
simultaneously jumps from "0" to "1," as shown in Figure 5(b). If the incident 
light (bias and signal) is sinusoidally modulated, i.e., the incident light 
power is varied from 0 to 12 mW sinusoidally, the bistable curve for every 
point of the 4 x 2 input array for OBD, is as shown in Figure 5(c). The photo- 
graph indicates that when the bias is 9 mW and the signal is 3 mW, all points 
in the 4 x 2 OBD, output array are at their high states and the output is "1." 
Using M, to block off S,S,, the 4 x 2 output array of OBD, is still at the high 
state, i.e., a bright spot "1" as shown in Figure 5(d). This illustrates that 
a 4 x 2 input array of "1" is written into and stored in memory. However, 
there is no light signal output from OBD,. When the optically stored informa- 
tion is needed, screen M, shown in Figure 3 is removed to allow the bias H, to 





(b) 
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shine on the 4 x 2 array in OBD, in the form of array P,P,. The stored data can 
then be read immediately because every point in OBD, jumps from a low state to 
a high state and the output optical signal is "1." 




















If we want to store the information from the __) Lise zero H, 
next interval in the same point of memory, f). 

the data stored must be cleared first. Before writing pulse 
storing new data, the bias light must be output of 0 
removed from the point to be cleared. After- — =. 
ward, bias light can be restored to begin [7 reading pulse 
storing new data. The writing, reading, and 

clearing of data are shown in Figure 6. | output of OBN, 
From this experiment, we know that the system Figure 6. The Process of 

can write, store, read, and clear data. In Writing Into and Reading 
optical computing, it can serve as a parallel From the Optical Cache 
cache memory array, as well as a parallel Memory 


random access memory. 
Conclusions 


The experimental optical parallel cache memory array consists of two uniform 
normal bistable ZnS devices. The optical system includes a high-efficiency 
planar thin-film beam splitter to provide a uniform beam array. A multipurpose 
SSR prism ensures the stability of the optical system. This simple optical 
system can perform basic parallel functions such as write, store, read, and 
clear in a random fashion. The simultaneous access of an optical signal array 
is primarily determined by the homogeneity of the beam array and the uniform 
ity of the bistable device. Various parameters of the optical cache memory 
array depend upon the upward or downward transition time. The reliabi:.ty of 
optical storage depends upon the stability of the bistable zone. In digital 
optical computing, this system can be used as a parallel cache memory array of 
a parallel random access memory array. 
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[Text] Abstract 


Matrix-vector and matrix-matrix multiplication have been achieved using four- 
wave mixing in LiNbO,:Fe crystal. Light intensities representing matrix 
elements are below 100 uW/mm? and the response time is of the order of seconds. 
The photorefractive effects affecting the accuracy of the operation are also 
discussed. 


I. Introduction 


Because all-optical computation based on photorefractive nonlinear optics is 
parallel and can handle data over a wide frequency band, it is a very attrac- 
tive subject. Yeh and Chiou’ first accomplished pixel-pixel multiplication 
using four-wave frequency mixing in BaTiO, in 1987. Then, these pixel products 
were summed with a cylindrical lens to complete the matrix-vector multiplica- 
tion experiment. This kind of vector-matrix multiplication can be extended to 
matrix-matrix multiplication by using composite light.* Matrix—matrix multipli- 
cation can also be done by finding the spatial convolution of four-wave 
frequency mixing. Chiou and Nguyen’ implemented this matrix-matrix multiplica- 
tion scheme with BaTi0,. As far as optical computing is concerned, the light 
intensity response of the photorefractive crystal must be sufficiently low so 
that a large number of light beams can be produced by a laser to complete the 
multiplication of higher-order matrices. Therefore, finding a crystal that has 
a low response to light intensity is critical. It was experimentally found 
that the optical response of LiNbO, with 0.1 wt (weight percent) Fe is below 
100 pW/mm?. Furthermore, the photorefractive effects, i.e., light-induced 
scattering’ and self-defocusing,* have substantial impact on accuracy. 
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II. Matrix-Vector Multiplication 
1. Multiplication by Four-Wave Mixing in Photorefractive Crystal 


The four-wave mixing process in a photorefractive crystal can be described by 
a steady-state wave-coupling equation. As shown in Figure 1, four beams of 
light with frequency w and identical polarization intersect in the 
photorefractive crystal. Let the photoelectric field of these four beams be: 








Ej(r,t) = Aj(r)exp[i(kj-r-wt) ] + c.c., (1) 

where j = 1, 2, 3, 4. Furthermore, the wave vectors A, z=l1 2 
satisfy these relationships: k, = —kz, k3; = —k,. a, As 
These four beams induce four refractivity gratings aN 
in the crystal. Two of them, i.e., transmission A 
grating and reflection grating, are the medium of 4 26 
four-wave mixing. The wave-coupling equations en) 
resulting are as follows®: 

= crystal” | AY 
(1) The four-wave mixing coupling equations axis 
primarily due to the transmission grating are: Figure 1. Diagram of the 


Four-Wave Mixing 











Sa - - LUA AL AIA) Ana Ay, | 
* - ae (A,A(+AzA,) As+aAz, 
ss = Ff (AAG +AGAy) Ayres, ' (2) 
a = DL (aaieara)ar-aal, 
_ 7 
scoot TY 





where n; is the refractivity grating amplitude, ¢, is the spatial phase 
difference between the grating and the interference fringe, c is the speed of 
light, and @ is the angle between the light velocities. 


(2) The four-wave-mixing coupling equations primarily due to the reflection 
grating are: 














dA eae 

=z * $ (nas +AzA,) A,-@A,, 

dA; a oi. . . * - 

dz = Tp (A,A; +A, A,) Ag +@A, ’ 

dA; . . ° ° 

Ge TT (MAS SAIAD AT oaas,| (3) 
dA eae 

= © ZT (nas +A, A,) A,-@A,, 

_ _ ion, exp (16,,) 
Y 2c cos 6 
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In LiNbO,:Fe crystal, the primary mechanism for electron transport is 
diffusion. Therefore, ¢, = ¢;; = (*/2).’ Because 7 is a real number and when 
there is almost no absorption (a = 0), these two sets of equations can be 
rigorously solved.® 


When equation (2), i.e., transmission grating, plays a dominant role, the 
conjugate reflectivity is 


| -2cT |? (4) 


R= A, (0) = 
Al(0)| | AT+(A2+4|c]?) 2/2] * 











A,(t) = 0, I,+I, = const = I,(0)+I,(0), I,+I, = const = I(t), 
A = (I,+I;) - (I,+I,) = I, (@) - I, (0) - I,(0), 
I,=I,+1,+1,+ T= 1,(0) + 7,(0 + 10), 
QO = (A*+4 |Cl|?)’/2, p = (yQ/2I,), T = tanhpe, 


' (5) 





|c|? is determined by the following equation: 
[JC]? — 1,(0)1,(2)] |AT + (A2+4]C]2)*/?7)2 + 4/C]7)/T)/71,(0)12(2) 
+ 2/C|21,(0) (A7+4|C|2)1/2(T+T") = 0. (6) 


Figure 2. Dependence of the Phase- s.b. experimental data 
Conjugate Reflectivity on 
the Pumping Intensities. The m 

‘ 

theoretical fitting curves S 2.54 
are obtained for y2£ =-1l, = 
I, = 1.0 mW. The numerical 
values above each curve cor- 
respond to the intensities 0 I (x100uW) 0.6 
of I, (x100 pW). e 











For I, = 1.0 mW, I,, I, < 100 pW, theoretical fitting curves for R versus I,(0) 
and I,(2) have been computed and are shown in Figure 2. When I, >> I,, Iz and 
I,, we have 


4 : -I,(0), Q . I,(0), a - (y /2),I, - I,(0), (7) 
Equation (6) can then be simplified as 
[Cj? = [(1-T)/(1+T) }#1,(0)12(2). (8) 


Furthermore, from equation (4), we get 





. |_1-exp(-y®) ; - 9 
R [EHO] 1,100 (0 I, (0) I,(0. (9) 


When equation (3), i.e., reflection grating, is dominant, its conjugate 
reflectivity is 
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ri), : 


tanh? [|g,|yp(é)] = 





R= 


I, (0) I,(0) [1-exp (-y@) }? (10) 


{I, (0) + [I (#) +I, (0) ] exp (-20) }#+Z, (0) Z,(0) [1-exp (-y@) ]? "| 








For LiNbO;:Fe, |y2| ~ 1. Thus, 


I, (0) [1-exp (-y@) }? 
I, (0) exp (-2y8) , (11) 


. |__1-exp(-y8 ‘ _ 
R T, (0) exp(-yo I, (0) 2, (0) « I, (0), (8). 





tanh? [|g,|n(@)] « 


~y 








4 


From equations (9) and (11), we can see that either transmission-grating— or 
reflection-grating—induced four-wave mixing can accomplish multiplication as 
long as the probe ratio q = {1,(0)[1,(0)+I,(2)]} >> 1. When 72 < 0, the 
transmission grating plays the major role. However, when yf > 0O, the 
reflection grating dominates. The use of yf is determined by the crystal axis 
and polarization of the coupled light (i.e., the effective photoelectric 
coefficient ye¢¢°®) . 


2. Matrix-Vector Multiplication by Image-Beam/Image-Beam Multiplication 


























photorefractive 
0, crystal uy 
¥ . 
) 1 
bes B \ » « ' 
TY 4 
BS , 2 
“21 
probe beas 
° (plane wave) 
U"U,U, im 
“4 | i coe Figure 4. Experimental Setup for 
Figure 3. Conceptual Diagram of Matrix-Vector Multiplication 
Matrix-Vector Multiplication M,, M2, M; are mirrors; L,, Ly, Ly, are 


cylindrical lenses; BS, BS,, BS, are 
beam splitters. CL is the collimator. 


Figure 3 shows the principle of matrix-vector multiplication. The matrix and 
vector elements are coded by two aperture masks with light passing in opposite 
directions. A collimated intense beam of light is used as the probe. Experi- 
mentally, a 1.5-mm-thick LiNbO,:Fe crystal is used for the image—beam/image- 
beam multiplication. These products are "read" by the probe light, which is 
focused by a cylindrical lens. The summation is done on the focal plane to 
complete the matrix-vector multiplication process. The experimental setup is 
shown in Figure 4. Figure 5 shows the result of the product of a 2 x 2 matrix 
and a 2 x 1 vector. For simplicity, a thick beam was used in the experiment, 
instead of using cylindrical lenses L, and L, to expand the bean. 
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Figure 6. Conceptual Diagram 
of Matrix Multiplication 
Using Convolution 








Figure 5. The Operation of 
Matrix-Vector Multiplication 


III. Matrix-Matrix Multiplication 
1. Convolution and Correlation Using Four-Wave Mixing 


White and Yariv proved that the correlation and spatial convolution of two 
patterns could be done by four-wave mixing.'® The crystal used in that 
experiment is Bi,,Si0,, and its principle is shown in Figure 6. After Fourier 
transformation, the object plane input amplitude U,(x,y), U2(x,y), and U3(x,y) 
become u,, uz, and uz (u; >> u,, uz) on the confocal plane. These three beams 
are coupled in the crystal to produce four-wave mixing. The conjugate of u;, 
u,, is converted to U,(x,y) by inverse Fourier transformation using lens Lj. 
It can be expressed as follows?®: 


U,(x,y) = pU3(x,y) & [U,(x,y) @ U2(x,y)], (12) 


where ~ is a transformation constant, * and e are correlation and convolution, 
respectively. 


When U, and U, are the matrix masks and U; is an aperture (6 function), we have 


C(x,y) = U,(x,y) @ U,(x,y¥) = JJ%s(&.n) 0, (x-8, y-m) dan 43) 
= > By yD114/8 (x-8is- Bir Yi Nyy) ’ 
ij, 1’! 
U(x, y) = WU; (x,y) WLU, (x, y) @U, (x,y) ] 
= vf [0, (En) Clx+8, yen) ded 
° (14) 


. ¥/ [8 (E.m) Clx+k, yen) dean 
= ¥C(x, y) , 
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where (&¥, nj,) is the coordinate of element y 





















































ij on matrix mask U,. Equation (14) shows May? "2 boo f¥ 
that after passing through lens L;, the b wg “al 
output amplitude of the conjugate wave ne (a) 4° — x (b) 
produced by four-wave mixing on the focal *21 *22 rt " 
plane is the convolution of U, and U2. Ys ~2a 0 2a 
Pibre)) 27B2y 82144 *BQ2b2, Baby 
If U, and Uz are the matrices shown in Figure CO [ansbo> Boabz0t8o1b49 a2. 
7(a) and (b), then their convolution is as dee fay pda, 844D14¢ yoboy aod, 
shown in Figure 7(c). The center column of bre | a,,bo5 8145, 2%8) >> Roltd 
the convolution corresponds to the elements (c) 
of UU. Figure 7. The Convolution of 
U, and Uz 


2. Matrix-Matrix Multiplication by Spatial 
Convolution 


Figure 8 shows the experimental setup. The horizontal distance of holes on U, 
and U, is 2a = 1.5 mm, and the vertical distances are 2b = 3 mm and 2c = 1.5 
mm. Figure 9(a)~(c) shows the experimental results of matrix multiplication. 
[The cross terms of the convolution (left and right column in Figure 7(c)) are 
not shown because they have no contribution to the matrix multiplication. ] 


01 21 
0 2) Gol =U 2] over 
U, Uy UU exposure 








Figure 8. Experimental Setup of Matrix- 
Matrix Multiplication 
M,, M2, M3; are mirrors; L,, L,, Ly are 
spherical lenses with equal focuses; 
U,, Uz, U, are matrix masks 





Figure 9. The Operation of 
Matrix—Matrix Multiplication 


IV. Experimental Results and Discussion 


Major factors affecting accuracy are self-—defocusing and light-induced 
scattering. Self—defocusing diffuses the output light, as shown in Figure 5. 
For long exposure, the light spots become undistinguishable. Self-defocusing 
is caused by an index change induced by a Gaussian beam, as shown in Figure 
10(b). If a light beam is uniform and weak, then the index of refraction 
change is flat and small. Consequently, the self—defocusing effect is negligi- 
ble. Light-induced scattering causes the conjugate wave to lose some energy. 
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As shown in Figure 10(a), multiplication is 10. 














accurate in the stable region marked in the D: beam 
figure. If scattering is strong enough, this | B Genstep 
stable region will disappear. Light-induced jp” 

scattering is dependent upon the diameter 2 

and intensity of the beam. Zhang Guangyin?} r lee D _ 
pointed out that there is a threshold dia- . Pee 2 net, cine 


meter for this type of scattering. There- (a) (bd) 

fore, a low-intensity and small-diameter 

beam can significantly reduce the effect of 

light-induced scattering and self-defocusing. Figure 10. Influences of Self- 


In this experiment, the image beam intensity Defocusing and Light-—Induced 
is less than 100 pW/mm* and the diameter is Scattering 

less than 0.1 mm. A stable area of more than (a) The energy loss caused by 
two seconds was obtained experimentally. In scattering. 

addition, in matrix-matrix multiplication by (b) The refraction index 
convolution, long-time exposure will lead to change caused by laser 
noise diffraction spots, as shown in Figure 9. beam. 1. Index change 
These noise diffraction spots are higher-order caused by Gaussian beam. 
diffraction points of the transmission or re- 2. Index change caused by 
flection grating. Hence, whether for matrix- weak uniform beam. 


vector or matrix-matrix multiplication, there 
is a definite time within which the result is 
accurate. 


V. Conclusions 


All-optical computation using photorefractive nonlinear optics has a wide 
range of future applications. In order to make it practical, the light 
intensity response of photorefractive materials must be significantly reduced. 
Experimentally, it was determined that LiNbO, crystal heavily doped with Fe 
(0.1 wt%) has a great deal of potential in this field. 


The authors wish to thank Associate Professor Zhang Chunping [1728 2504 1627} 
and Instructors Li Laobing [2621 5071 0365] and Cui Yufen [1805 3768 5358] for 
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[Text] Abstract 


An IPA (Inter-Pattern Association) memory based on Fresnel hologram is 
presented. It has higher recognition capability for nonindependent storage 
modes. Using Fresnel holograms, an optical IPA memory with large number of 
neural elements can be implemented. 


I. Introduction 


Neural networks have been an active area of research in recent years. Associa-— 
tive memory is the most critical component in a neural network. In most 
optical neural networks, Hopfield’s model is used to form associative memory. 
However, Hopfield’s model does not consider the relation between various modes 
in storage. It only stresses the inter-pattern association of elements. 
Hopfield’s model is valid when storage modes are independent and the ON and 
OFF probability for each element is equal. Such a system is a poor model for 
dependent modes. When it is necessary to distinguish similar modes, the 
difference between them becomes very important. Therefore, when we develop an 
IPA memory, the similarities and dissimilarities of storage modes must be 
considered to make them function differently in the association process. 


In 1990, Taiwei Lu’ proposed using an IPA model to fabricate an associative 
memory. An optical neural network, consisting of a monitor, lens array, spa- 
tial light modulator, and computer, has been constructed. The IPA associative 
memory matrix and the effect of the lens array on each pixel are shown on the 
monitor. The results indicate that the address—finding capability of the IPA 
neural network is substantially higher than that of the Hopfield’s model. 
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Especially when the Hamming distance of the storage modes is small, it is 
capable of finding the address correctly. Thus, the IPA neural network has 
more stable capability than that of the Hopfield’s model. However, because of 
the resolution of the monitor, the number of pixels stored in the system 
cannot be very large. 


A hologram associative memory is particularly suitable for the storage of two- 
dimensional pictures with multiple pixels. A theoretical analysis of the 
characteristics of Fresnel holograms is presented in this paper. An IPA 
associative memory was fabricated on the basis of this analysis. Furthermore, 
a large-pixel IPA optical associative memory has to be implemented. 


II. Theoretical Analysis 


The IPA model decomposes every mode into several parts based on how one mode 
coincides with other modes in the set of storage modes. Each part has a 
different function in the associative process. A unique part of the mode has 
a more important function than other common parts in association. 


For example, in storing three modes, when some part of the input address- 
finding mode is ON, if the corresponding parts of the three modes are all ON, 
then the input mode may be any one of the three modes. There is no way to tell 
which one; we can only say there is an input mode. Hence, this part only 
stimulates the common parts. If the corresponding parts of two modes are ON 
and that of the third mode is OFF, then the input mode is one of the former 
two and it is not the latter. However, it is not possible to distinguish 
between the two. Then, this part stimulates the common part of the former two, 
but suppresses a unique part of the latter. If only the corresponding part of 
one mode is ON and those of the other two are OFF, then the input must be the 
former. This part stimulates all the parts of this mode and suppresses some 
part of the other two modes. Taiwei Lu proposed a triple gray-scale associa- 
tive matrix. The elements are -1, 0, and 1. In an optical neural network, it 
is easy to implement stimulation and suppression. This kind of associative 
matrix is relatively easy to materialize optically. Nevertheless, the function 
of each part in the associative process in an IPA memory is proportional to 
the number of pixels contained in this part. When the Hamming distance of two 
storage modes is small, there is a substantial overlap between modes. Their 
unique parts play a limited role in the associative process. It is more 
difficult to resolve these modes. We believe that we should add a gray scale 
in the associative matrix to balance different parts. Since the system chooses 
a threshold for output, the gray-scale requirement is not very stringent, 
which is easier to realize optically. 


The following is a theoretical analysis of IPA memory with Fresnel hologram. 
The optical path used to record the Fresnel hologram is shown in Figure 1. P, 
is the imaging plane of plane P, through lens L,. P; is the imaging plane of 
plane P, through lens L,. During recording, a transparent plate with a trans- 
missivity of a(p,q) is placed on plane P, to come in contact with the holo- 
graphic mask. It is illuminated by another focused light. Another transparent 
plate with a transmissivity of b(x,y) is placed on plane P,. A diffuse scatter- 
ing screen is placed in front of it. 
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Let us assume that the diffraction field of b(x,y) on plane P, is B(p,q). In 
addition, the distribution of this field on plane P, is noted to be D(x’, y’) 
b//2(x', y'). Then, we have 


B(p,q) = [ [oily bi? (x,y) exp [-ik(d+p sin @) ) (1) 
‘exp [ik{(p cos 6-x’)?+(q-y’)?} /2d] dx’fy’. 
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Figure 1. Fabrication of Memory — s 
Matrix for an IPA Holo- Fomput er 
graphic Associative Memory Figure 2. IPA Holographic 


Associative Memory 
The exposure on the recording medium is 


proportional to 


| (exp[ik(p?+q?)/2f] + B(p-q))a(p,q) |”. (2) 


Here, d is the distance between the holographic mask and plane P3, K is the 
wave number of the laser, f is the focal length of the lens, 6 is the angle 
between the optical axes of the two lenses, and D(x’, y’) is the diffuse 
scattering screen. Masks representing modes aj and their corresponding 
stimulating and suppressing parts b; are placed in positions a and b, 
respectively, prior to exposure. After processing, the amplitude of the 
transmissivity of the hologram is proportional to the amount of exposure. 
Terms contributing to first-order holographic diffraction include: 


t(p,g) = > {exp [-ik(p?+q?) /2f)] B,(p, @) }a, (p,q) . (3) 


Figure 2 shows the optical path of a holographic IPA memory. After placing the 
recorded Fresnel hologram on plane Pz, when the input object mask a’(x,y) is 
placed on plane P, and is illuminated by a collimated beam, the transmission 
holographic optical field is 


exp[-ik(p*+q?)/2f]a’(p,q)-t(p.q). (4) 


On output plane P;, the optical field can be expressed as: 


y [fase aa'ip, q) B, (p, gq) exp [-ik(d+p sin 6) } (5) 
‘exp[ik{(p cos 6-x’)? + (q-y’)?}/2d): dpdq. 


31 








Let Fj(x,y) be the Fourier transform of aj;(p,q)-a'(p,q). Then, the output 
optical field can be written as: 


> flr (x-x!, y-y’) D, (x’, y!) bi? (x!, y’) 
exp [ik(x?-x+y?-y”) /2d] > dx/dy’. 


(6) 


Because of the presence of the diffuse scattering screen Dj(x', y’), the 
intensity distribution can be expressed as: 


yf ftaias - 8(x-x!, y-y') by (x!, y’) dx'dy’ 


(7) 
= >> {a,a'}b, (x,y) . 


The bracket ({aja’) represents the inner product of g and a’. Here, the diffused 
light spot on the output plane is neglected. In the associative process, inner 
products of the input object and different parts are obtained to determine the 
corresponding stimulation and suppression distribution. The weighted associa- 
tion is the magnitude of the inner product. The magnitude of the associative 
matrix elements can be implemented by adjusting the diffraction efficiency of 
the holographic mask. In this work, different diffraction efficiencies were 
obtained by adjusting the exposure time to record the holographic mask. 


III. Optical Implementation 


Based on the above theoretical analysis, an IPA Fresnel holographic memory has 
been realized optically. Experimentally, three 24 x 24 pixel English letters, 
X, Y, and Z, were chosen as the patterns to be stored in binary form. Based on 
this concept, pixels can be divided into seven categories: (XAY)UZ, (YNZ)UX, 
(ZAK)UY, (XNY)NZ, (YNZ)NX, (ZNX)NY, (XNYNZ)Ng, where N, U, and — represent 
logic operations AND, OR, and NOT, respectively. ¢ represents the set of 
pixels none of the three patterns uses. The matrix elements with suppressing 
effect are negative in value. This work also adopts a commonly used scheme to 
provide a bias for every matrix element so that the minimum occurs at zero 
transmissivity in order to compute negative values optically. 


In this experiment, a 30 mW He-Ne laser was used as the light source. An HP- 
366-1 holographic mask was used to record multiple exposures. Since there is 
no overlap between the light-transmitting pixels in every category, the 
maximum number of exposures at any point of the holographic mask is limited to 
one to ensure the diffraction efficiency of the mask. A two-dimensional CCD 
detector array is used to receive the output picture. The signal is input into 
a computer and the computer selects a threshold for the output picture. The 
associative result is then shown on the monitor. Figure 3(a) shows the picture 
stored in the system. Figure 3(b) is a part of the input stored by content- 
address finding. Figure 3(c) is the associative result with this kind of 
input. An IPA memory with Fresnel hologram can be used to implement pattern 
association with a large number of pixels. Experimentally, we have demon- 
strated it with 24 x 24-pixel patterns; however, this by no means is the 
limit. Because the resolution of the holographic mask is very high, it does 
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not limit the pixel number. The maximum number of pixels is determined by the 
input mask. 





(a) (b) 





Figure 3. Stored Patterns and Recalled Outputs 


Conclusions: An IPA memory with Fresnel hologram has been fabricated to 
increase the number of pixels to be stored in a pattern. In addition, the gray 
scale in the associative matrix is increased in the IPA memory to ensure high 
resolution of storage modes with small Hamming distance. 
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[Text] Abstract 


A hybrid optical—digital system capable of solving a set of linear equations 
using an iterative method is presented. The optical part primarily consists of 
a single holographic mask used to carry out vector-matrix multiplication. The 
remaining part, including a CCD detector and a microcomputer, is for vector 
measurement and addition. This hybrid system was used to solve a set of four- 
dimensional equations. Compared to theoretical solution, the error of the 
experimental result is approximately 5 percent. 


I. Introduction 


In recent years, there are a number of studies’° on optical computation of 
matrix multiplication. However, optical solution of linear equations is still 
rare. In 1982, Casasent presented a scheme for solving linear equations by 
using an optical vector-matrix multiplier with acousto-optic modulation and 
electronic feedback iteration.? In 1986, Caulfield, et al., presented a dual- 
mode optical scheme for solving linear equations using an analog and a digital 
optical computing system.‘ Nevertheless, they failed to present any experi- 
mental results in their papers. 


It is well known that there are numerous iterative methods for solving linear 
equations. The most common is a simple iterative scheme. This scheme involves 
the inner product of a vector and a matrix and a vector addition in every 
iteration. The former can be done optically. As for the latter, it is more 
appropriate to use an electronic system. Therefore, using a hybrid optical- 
digital system to solve linear equations iteratively is a very logical choice. 
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II. Principle 
Let a set of linear equations be 
AX = B (1) 


where A and B are a two-dimensional matrix and a one-dimensional vector, 
respectively, and X is an unknown vector. After a primary transform of A, the 
iterative equation for this set of equations is 


x (k+l) _ Tx‘*) + D (2) 


where matrix T is related to A, and vector D is a function of B. The 
superscripts k and k+l represent the number of iteration. Equation (2) 
converges only when one of the following conditions is met: 


max d le(4.5)1 < i, | 
i¢gign el 


a 
max > lei. 1 <1, YO leyl <1, 
=1 i,j*l 


1¢jén 


(3) 





where t(i,j) is the element of matrix T in the ith row and jth column. If none 
of these three conditions can be satisfied, it is necessary to transform (2) 
further to alter matrix T and vector B until t(i,j) meets one of these 
conditions. 


The accuracy of solving linear equations using (2) can be estimated as 
follows: 


xckrt) a x(k | < €, (4) 


where ¢« is an arbitrary small quantity. 


holographic 
The hybrid optical-digital system designed — (3 
and constructed in this work is shown in lens - o 
Figure 1. In this system, the complex ~ 7 tf 
amplitude distribution of the modulating 
holographic mask H(m) is defined by the ‘ ‘ 1 . . 











following equations®®: x, x 


2 
input plane ant plem output plane 
Figure 1. Scheme of Optical- 


t(i,j) = > Gi, m) H(m) G(m, j), Digital Computing System 
1 = is 2, we, N,, me i, 2,°, N, (5) 
j = 1, Zs. ee N, 


In these equations, G is the Fourier transform factor; N,, N, and No are the 
number of sampling elements on the input plane, mask plane, and output plane, 
respectively. When the following equation is satisfied by optical parameters 


(AxAx;/fA) = N*, (AxAx2/fA) = N27, N = N,No, (6) 
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equation (5) has a unique solution. Moreover, based on unitarity, it is 
possible to obtain the following solution directly: 


y &, 


H(m) = ) > 2, O° lm 1) (1, JG° (Jom /sin c(r,m/N), (7) 
=1 #1 





In equations (6) and (7) , Ax,;, Ax, and Ax2 are the sampling periods of the 
input, holographic mask, and output phase, respectively, f is the focal length 
of the Fourier lens, A is the wavelength, and r, is the ratio of input plane 
sampling element aperture to sampling period. Therefore, with a given matrix 
(from iteration equation (2)), the complex amplitude distribution of the 
holographic mask H(m) can be determined in equation (7). The design parameters 
of the optical system can be selected over a reasonable range based on 
equation (6). 


III. Experiment and Results 


The following is an arbitrary four-element set of linear equations: 


38x, + 3X, + SX, + 14x, = 70, 100xX, + 91x, + 26x, + 48x, = 300, (8) 
9x, + 6X, + 31x, + 7X, = 115, X, + 17X, + 30x, + 94x, = 200, 


This set of equations (8) is transformed into an iterative format (9): 














(x,)"°") /-0.24 0.06 0.10 0.28\(%)" (1.4) 
X, __| 0-24 -0.15 0.16 0.20]/%] | /1.6) (9) 
x, 0.18 0.12 -0.38 0.14] |x, 2.3 
x] ( 0.01 0.17 0.30 -0.06, X,) (2. 0) 














Comparing to equation (2), we get: 


(-0.24 0.06 0.10 0.28) 
Ts 0.24 -0.15 0.16 0.20 (10) 
0.18 0.12 -0.38 0.14 


. 0.01 0.17 0.30 -0.06, 








D’ = (1.4, 1.6, 2.3, 2.0)'. (11) 


(T represents transformation.) After verification, matrix T meets the conver- 
gence equation (3). Equation(10) shows that N, = N, = 4 and N = 16. Substitute 
matrix T into (7) and let r, = 0.4, then we get the complex amplitude distribu- 
tion H(m) of the holographic mask. Its numerical values are shown in Table 1. 


The design parameters of the optical system are determined based on equation 
(6) and the existing experimental conditions. They are Ax, = 0.2052 mm, Ax, = 
0.8206 mm, Ax = 0.0800 mm, f = 415 mm, and A = 6.328 x 10°‘m. 


Figure 1 shows the experimental scheme of the hybrid optical-digital system. 
The holographic mask on the confocal plane was developed using computer 
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generated holograms based on the data shown in Table 1. The amplitude and 
phase distribution extend along the x axis. Each amplitude phase value is 
represented by a rectangular grating extending along the y direction. If the 
grating period is p, starting position is y and the light transmitting aper- 
ture is Ay, then the first-order diffraction phase of the grating is ¢ = (y/p) 
and amplitude A = [sin(xAy/p)/(x/p)]. In this experiment, the grating period 
is 0.1 mm. The positioning accuracy of the grating pattern generator (GCA3600) 
is 0.1 wm. Therefore, the phase and amplitude of the holographic mask can be 
kept at one ten-thousandth.’ The magnitude of the input vector is represented 
by the sampling element aperture size on the input plane. The output vector is 
measured by the CCD detector array. The result undergoes an A/D conversion and 
then is input into a microcomputer. Vector addition is done in the 
microcomputer. 


Table 1. Amplitude—Phase Values of the Holographic Mask H(m) 












































1 a 3 4 5 * 7 . 
Amplitude 6.211 1.008 0.203 | 0.199 | 0.704 | 0.543 0.657 0.622 
Phase 1.222 | -3.057 | -1.923 | 0.828 | —0.193 — 0,300 | 2.033 — 0.088 

” 9 10 | 1 | 12 | 13 | 14 | 15 16 
Amplitude 0.622 | 0.657 | 0.542 | 0.704 | 0.192 | 0.203 | 1.000 0.311 
Phase 0.083 | -2.033| 0.300 | 0.193 | -0.8a2| 1.993 | 3.087 ~ 1.222 











The position of the holographic mask is adjusted in the following manner. On 
the z axis, the confocal plane is found by using a ground glass to find the 
maximum spot. A center line is set along the x-axis. It is lined up with the 
center line of sin c diffraction generated on the holographic mask plane by 
the input aperture. Alignment along the y-axis is not important because the 
grating is periodic. 


Let us assume that (1 1 1 1)' is the initial solution xX") of the set of 
equations. An input plate representing X‘°’ is placed on plane x, shown in 
Figure 1 and the result obtained by the optical-digital system is X“’’. xX‘? is 
then used as the input vector for the hybrid system in the second iteration. 
After repeated iterations, the solution of the set of linear equations can be 
obtained. Table 2 shows the experimental data obtained by the hybrid system 
and the theoretical values derived from computer-simulated iterations. They 
are in very good agreement. 


The real solution solved by a computer is X" = (1.01, 0.79, 3.04, 1.01)'. After 
six iterations, the solution obtained by the hybrid optical-digital system has 
an error of 





2 
x(n) -x(n) 
x’ (n) 


2 5.2%. (12) 
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Table 2. The Iterative Solutions of the Four—Dimensional Equations Produced by 
the Hybrid Optical-Digital System and Theoretical Calculation 






































solutions (1) solutions (2) solutions (3) solutions (4) 
iterative time : = : : - _ : “a 
exp. theo. exp. theo. exp. | theo. csp. theo 
1 1.20 | 120 | 103 | 115 | 296 | 2a¢ | 149 | 158 
2 1.04 0.95 0.81 | 0.81 2.44 | 2.57 120 | 1a 
3 0.99 | 09s | 081 | 0.88 | 2.78 | 2.8 | 1.95 | 135 
‘ 0.92 0.97 0.77 j 0.81 2.95 2.94 1.11 | 1.06 
5 1.03 | 0.99 | 0.77 | 0.80 | 2.95 | 2.90 | 114 | 1.03 
6 1.05 1.00 0.82 | 0.80 8.02 | 3.02 1.10 | 1.08 











Due to the parallel nature of the optical system, linear equations can be 
quickly solved using a hybrid optical-digital system. However, since it 
involves an analog operation, its accuracy is somewhat limited. Nevertheless, 
it is useful as the front end of the scheme proposed by Caulfield for using a 
dual-mode optical computer to solve linear equations, because it only requires 
a rough solution at high speed. 


The authors wish to thank Zheng Shihai [6774 1597 3189] and Li Dehua [2621 
1795 5478] for their assistance. 
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[Text] Abstract 


For the first time, multiple-valued optical logic has been demonstrated using 
a bistable laser diode (BILD) or bistable light emitting diode (BILED). Basic 
ternary logic functions such as COMPLEMENT, MINIMUM, MAXIMUM, CYCLE, and 
LITERALS are experimentally demonstrated. A triplexer circuit design is 
presented. 


Hybrid optoelectronic logic circuits have attracted a great deal of attention 
lately.1-3 A BILED-based hybrid optical multiple-valued logic circuit is 
presented in this paper. Experimentally, five basic ternary logic functions in 
POST algebra, i.e., COMPLEMENT, MINIMUM, MAXIMUM, CYCLE, and LITERALS, have 
been implemented. Furthermore, a ternary T-gate, or triplexer, design is 
shown. 


Compared to binary logic, a multiple-valued logic variable carries more 
information. Therefore, logic operations carried out by a multiple-valued 
logic gate require less logic gates and connections.‘ Just as Boolean algebra 
is the theoretical basis for binary logic, there is a theoretical framework 
for multiple-valued logic. POST algebra is a common one. Considering ternary 
logic where R = 3, its logic variables are xe{0,1,2}. The basic logic 
functions in ternary POST algebra include: 


COMPLEMENT xX=2=2-x> 
MINIMUM Min (x,y) = x, if x < y. 
MAXIMUM Max (x,y) = x, if x2 y. 
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CYCLE Suc(x) = x + 1, (mod 3) 


nb_{ 2 if asxsb 

LITERALS “x { 6 elsewhere 

The multiple-choice switch T(x,y,Z:K) = : if ret 
z if K=2 


or T-gate operation is: 


where COMPLEMENT, MINIMUM, and MAXIMUM are essentially equivalent to NOT, AND, 
and OR in binary logic. CYCLE can be used to implement the transformation of 
a logic variable. LITERALS represents the state of a logic variable. If a = b 
= K, then LITERALS provides the state of x at 0, 1, 2 as Jo(x), J,(x), and 
J2(x). The T-gate can be used to select the output of a variable K at different 
values. 

+V 
There are two factors to be considered 
in implementing multiple-valued logic. Re 
(1) Are there multiple input elements 
that can provide multiple output? 
(2) Can elements’ be_ conveniently 
connected? It is simple to connect ele- 
ments in an optoelectronic circuit. We 
used a tristable circuit based on two 
BILEDs, as shown in Figure l(a), to 
create a unit element for the optoelec-— 
tronic ternary logic. In the figure, the 
LED is a red GaAsP light emitting diode, 
PD is a 3DU32 silicon phototriode, and T 
is a transistor; f, and fz, are the two 


optical feedback coefficients for the Tristable Circuit or a Photo-Current 


two bistable circuits. The output ampli- Mirror; (b) Output-Input Character- 
tude is controlled by bias V and resis— j¢ric Loop 


tances R, and Rj (i = 1,2). The thresh- 

olds to turn on BILED, and BILED2 are 

controlled by resistors Rg, and Rez. If Re, > Rez, then Pon; < Pon2, where 

Poni and Pon, are the thresholds to turn on BILED, and BILED2, respectively. Due 
to the branching of the circuit, a tristable output is achieved with two 
saturated currents and the switching of BILED, and BILED, at different 
thresholds, as shown in Figure 1(b) [photograph not reproduced]. 








Figure 1. (a) Schematic Diagram of a 


If the three steady states are 0, 1, and 2, and the switching thresholds 
satisfy the following: 


l < Poy < 2 (1) 


then the circuit can have a ternary output with respect to a ternary input in 
a reproducible manner. This is analogous to the current mirror in the I?L 
multiple-valued logic design in electronics. The tristable circuit shown in 





Figure l(a) is called a light-current mirror, a device which plays a key role 
in almost all ternary logic circuits. The following is a description of 
optical ternary logic gates and some experimental results. 


Figure 2(a) shows the circuit of a ter- 
nary COMPLEMENT gate, which consists of 
a light-current mirror and a light 


emitting diode LED,. Because the output | Ry 
of LED, is opposite to that of the 
light-current mirror, its output P, — j— f= } — a | 
appears in a downward tristable form AaB, 
with a triangular wave input P;. When BILED, BILED, , 
x = 0, the light-current mirror cannot 
be turned on and all the current flows = + 4 

(a) 


through LED,; consequently, Pp) = 2. When 
x = 1, only a part of the current flows 
through the light-current mirror, and 
output P, = 1. When x = 2, all the 
current flows through the light-current 
mirror and P, = 0. The input and output 
waveforms are shown in Figure 2(b) [photograph not reproduced]. 


1. COMPLEMENT Gate i V 

















Figure 2. (a) Schematic Diagram of a 
Ternary COMPLEMENT Gate; (b) Input 
and Output Waveforms 


When the amplitudes are 1 and 2, the input and output current are expressed by 
the LED current, which is 9.5 mA and 19.0 mA, respectively. When Py, = 0, 
because of the added clamping diode D, the output of LED, is essentially held 
to 0. At this time, the current passing through LED, is approximately 0.8 mA. 


2. MINIMUM Gate 



































+V 
A ternary MINIMUM gate consists of two light- R, 
current mirrors, as shown in Figure 3(a), LED, —= 2 
where BILED, and BILED,, and BILED, and BILED,, ‘om + 
form two separate light-current’ mirrors. z a 
Furthermore, BILED, and BILED,, as well as — |BILED, BILED3 —— 
BILED, and BILED,, are also connected in ba 7 
series, respectively. The switching thresholds y ——~+ 4 bd 
of the bistable circuits are regulated by =~ |BILED, BILED, dn 
resistors Rg; (i = 1, 2, 3, 4): — 
t =z: = 
0 < Pon = Pon < 1 ) 
1 < Pons = Pong < 2 (2) 


Figure 3. (a) Block Diagram of 
a Ternary MINIMUM Gate; 


At this point, BILED, and BILED,, and BILED, (b) Experimental Results 


and BILED, are two separate elements to be 
controlled by x and y, respectively. When the 
input signals are 1 and 2, they are switched on simultaneously. If x and y are 
1, then BILED, and BILED, are switched on at the same time and the output is 
P) = 1. If x and y are 2, then BILED, and BILED, are also switched on and Py = 
2. Thus, it accomplishes the MINIMUM operation for x and y. Figure 3(b) 
[photo-graph not reproduced] shows the experimental results. 
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Figure 5. (a) Block Diagram of a 
Figure 4. Circuit Diagram of the Ternary MAXIMUM Gate; 
Improved Ternary MINIMUM Gate (b) Experimental Results 


In reality, a ternary MINIMUM gate can be accomplished with an improved light- 
current mirror to further simplify its structure, as shown in Figure 4. Two 
phototriodes PD, and PD, are added in series to the input end of two BILEDs. 
Since the photocurrent of the phototriode saturates at the minimum input light 
intensity, a single light-current mirror can act as a ternary MINIMUM gate as 
long as condition (1) is met. 


3. MAXIMUM Gate 


Figure 5(a) shows a schematic of a ternary MAXIMUM gate comprised of two 
parallel-connected light-current mirrors. The light-current mirrors are 
similar in construction to those described above. In addition, LED, and LED), 
and LED, and LED,, are connected in parallel. Therefore, BILED, and BILED,, and 
BILED, and BILED, form two separate elements to complete the MAXIMUM operation 
with amplitude greater than or equal to 1 and that for amplitude equal to 2, 
respectively. Once BILED, is switched on, BILED, is also switched on due to 
light feedback. BILED, and BILED, behave similarly to BILED, and BILED,. If the 
two light-current mirrors satisfy condition (2), then only BILED, and BILED, 
are on when x = 1 or y = 1, and Pp = 1. When x or y = 2, all the BILEDs are on, 
and P, = 2. Figure 5(b) [photograph not reproduced] shows the experimental 
results. 


4. CYCLE Gate 


Figure 6(a) is a schematic of a one-step-right CYCLE gate. A light-current 
mirror, comprised of BILED, and BILED,, is connected to BILED, in parallel. A 
ternary signal and a constant light with amplitude 1 are simultaneously input 
into the system. The switching threshold of BILED, meets the following 
condition: 


2 << Pons < 3 (3) 


Therefore, after the photocurrent jumps up to 2 twice in a row, BILED, is 
switched on to allow the current to flow from the light-current mirror through 
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very low relative to that of the light- 
current mirror). Figure 6(b) [photograph 
not reproduced] shows the experimental 


BILED, (if the resistance of BILED, is hee 





i 
2 
input signal is x = (210) and the output LED.Y ———P 









































results of the one-step-right CYCLE. The 
is Pp = (021). 
z xz Q 

Signal x can also be cycled one step ~~ | gyep » | eep, | 
left, 0 + 2, 17+ 0, 2 + 1. Based on L — | ao| SLED. 

, : | 
definition, one step left is equivalent L T } 
to two steps right. Therefore, one-step- = ee OU ¥ 
left may be accomplished by going _ ©) 
through the gate shown in Figure 6(a) Figure 6. Schematic Diagram of a 
twice. Nevertheless, it may also be done One-Step-Right CYCLE Gate 


with a one-step-left CYCLE gate. Figure 

7(a) shows a sketch of the gate. By 

adding BILED, to the circuit and changing the constant light amplitude to 2, 
the threshold for switching BILED, on becomes: 


3 << Pow < 4 (4) 
Hence, after BILED, is switched on, a portion of the current flows through the 


light-current mirror to result in P,) = 1. This operation makes it possible to 
cycle x one step left. 
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Figure 7. (a) Block Diagram of a One-Step-Left CYCLE Gate; 
(b) Circuit Design of a Ternary Half—Adder 


It should also be pointed out that if the constant light 2 is replaced by 
another signal y, the circuit shown in Figure 7 is also an adder without 
carry. When x+y < 2, the light-current mirror output is the direct addition 
result. When x = 1, y = 2, or x = 2, y = 1, BILED, is turned on and Py, = 0. 
When x = y = 2, BILED, is switched on and Pp) = 1. Hence, the operation x AND 
y is realized. Please note that when BILED, is switched on, the output of 
BILED, is not reduced to 0. Hence, the output of BILED,, P)3;, can be used as an 
input for another BILED where the threshold is very low and the output 
amplitude is 1. This circuit can be used to output a carry. Thus, we have a 
ternary half adder, as shown in Figure 7(b). 
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Figure 8. Block Diagram of 
LITERALS Gate 














5. LITERALS Gate 





Functions such as Jo(x), J,(x), and J2(x) can 
be generated by BILED circuits comprised of a 
binary NOT gate and an exclusive OR gate in 
parallel, as shown in Figure 8. BILED, and LED, 
form a binary NOT gate,°* while BILED, and 
BILED, form an exclusive OR gate. The output of LED,, BILED,, and BILED, are 
Jo(x), J,(x), and J2(x), respectively. When x = 0, the BILED circuit is not 
switched on and all the current goes through LED,, i.e., Jo(x) = 2 and J,(x) 
= J,(x) = 0. When x = 1, both BILED, and BILED, are turned on and current flows 
through BILED, and BILED,. Since there is a clamping diode in the LED, branch, 
the output of LED, drops to 0. In this case, J,(x) = 2 and Jo(x) = J2(x) = 0. 
When x = 2, the BILED is also switched on. Because the resistance of BILED, is 
far less than that of BILED,, and there is a clamping diode in the BILED, 
circuit as well, J,(x) falls to 0, i.e., Jo(x) = 2, Jo(x) = J,(x) = 0. 
Experimentally, the output of BILED, is far greater than 2. It is possible to 
normalize it to 2 by using a shunt resistor in parallel with LED,. 





Figure 9. Block Circuit of a 
Multiplexer or T-Gate 


6. T-Gate 
Let us re-express the T-gate as follows: 
T(x,y,z; K) = Min(J9(K), x) + Min(J,(K), y) + Min(J2(K), z) (5) 


Therefore, a ternary T-gate consists of a LITERALS gate and three MINIMUM 
gates. Figure 9 shows the block diagram of a T-gate. Through the LITERALS 
gate, a control variable K results in three functions J)(K), J,(K), and J2(K). 
These three functions and the three signals x, y, and z pass through three 
MINIMUM gates, respectively. Different variable output is obtained based on K. 
The simplified MINIMUM gate circuit shown in Figure 4 may be used as the 
ternary MINIMUM gate. 
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[Text] Abstract 


A programmable real-time optoelectronic hybrid threshold logic processor is 
presented. Its applications in optical computing and optical information 
processing are discussed. 


I. Introduction 


In optical computing, the most intriguing subject is the optical neural net- 
work, which mimics the neural logic processes of the human brain.’ Because the 
human brain can process logic in a highly parallel, fuzzy, multiple-valued, 
weighted, threshold-taking and feedback-controlled manner, it is necessary to 
investigate an optical logic stronger than binary vertex primitive logic, 
such as fuzzy logic, multiple-valued logic, and binary and multiple—valued 
threshold logic. In recent years, optical multiple-valued logic operations 
have been implemented in several optical systems.** Although Arrathoon? 
suggested using waveguides to implement optical threshold logic, very few 
simple optical threshold logic schemes are available. For the first time, an 
LED (light emitting diode) and PD (phototriode)-based optoelectronic hybrid 
threshold logic processor is presented in this paper. The processor can carry 
out real-time programmable control by optical signal to obtain two weighted 
and threshold values, i.e., both positive and negative weighted and threshold 
values. Therefore, the logic processor itself is a complete logic set. This 
logic processor is simple in structure, reliable, and more importantly com— 
patible with VLSI (very large-scale integration) technology. Optoelectronic— 
microelectronic integration can be directly implemented on a GaAs or GaAs/Si 
wafer. Hence, it has a bright future in optical computing and digital optical 
information processing. 


46 








II. Optoelectronic Hybrid Threshold Logic 
Processor 


Figure 1 shows the schematic diagram of a 
threshold logic gate. It has n binary input 


variables x;(i=-l,2,------ ,n). Each variable 
is accompanied by a weight factor 
wj(i|el,2,------ ,n). There is a _ relation 


between output y and expected threshold T®: 


{ a 
nh, for >> w,X,2T, a 








a 
0, for p> w,X,<T. 
“1 Figure 1. Schematic Diagram of 


Only when the weighted input variable is Threshold Gate 


greater than or equal to the threshold is 

the output y equal to 1. Otherwise, it is 0. If a switching function f(x) is 
a threshold logic function, i.e., it can be produced by a threshold logic 
gate, then there is a weight-threshold set (w,,w2,---,W,;T) that satisfies 
equation (1). Mathematically, if (w;,T) covers all real numbers, then a single 
threshold logic gate forms a complete logic set. This means that all logic 
functions can be implemented by a finite number of threshold logic gates. 
Furthermore, a threshold logic gate has more input capacity than a binary 
basic logic gate. Hence the number of logic elements and transmission lines 
can be significantly reduced. 
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Figure 2. (a) Circuit of Optoelectronics Threshold Logic 
Processor; (b) Part of the Circuit of Figure 2 
(a) Which Only Consists of Input x; 


Figure 2(a) shows the circuit of an optoelectronic hybrid threshold logic 
processor. In reality, it is a variation of a bistable LED circuit.‘ Experi- 
mentally, optical feedback was removed to increase the response speed of the 
logic gate. The LED output is the logic gate output y. Each weighted input end 
contains four silicon PDs in series. Two serve as the input of the input 
variable and the other two are the input of the weight signal. Each weight w; 
is divided into two parts: w;* and w; . When x; is positively weighted, 
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wj =~ wij’ and w; = 0. If x; is negatively weighted, then w; = -w; and w;* = 0. 
Threshold T is also controlled by two PDs connected in series. Similarly, it 
has a positive and negative value. Different weight and threshold values can 
be obtained by varying the intensities of the weight and threshold signal. 
Consequently, this threshold logic processor becomes real-time programmable. 


In order to illustrate how real-time weight and threshold can be implemented 
by the processor, let us use an example where there is only one input variable 
x;, as shown in Figure 2(b). The input x; illuminates PD;* and PD; and the 


weight control signals w;* and w; shine on pp*, and pp:,, respectively. The 
current intensity Igj can be expressed as: 


To, = Tas - Tas (2) 


Assuming all PDs are identical and operate in the linear amplification area, 
then the output photocurrent Ig is proportional to the input light intensity 


Pi, i.e., Igj = KPj;, where K is an optoelectronic conversion constant. Because pp} 
and pp;, (PDj and pp;,) are -onnected in series, 7}, (73,) is proportional to 
either the x; or w;* (w; ), whichever has the lower optical intensity; i.e.: 


Ij, = K min(P(x,), P(w{)) (3) 


where min(x) is the minimum threshold function. P(xj) and p(w?) are the 


intensities of xj and w}, respectively. The way weight and threshold values are 
changed is different from multiplying the amplitude of the logic variable by 
a factor wi. From equation (3) we can see that p(w?) is far less than the 


intensity at xj = 1. Hence, 73, is a function of p(wf) only when xj = 1. We 
have: 


¢ _ JKP(wj), for x; = 1 A 
Ts = 15 for x, = 0 (4) 


Since P(X4) |p0r>P (2) Iyer let P(X 4) |yi0r = gP(w}) luges® this means that the 
intensity amplitude at x; = 1 is evenly divided into q parts. If the weight of 
x; is tr(l < r sq), then P(wf) =r P(wf) luges 
50 wW and at wi = 1, P(wj) = 5 wW; then q = 10. If the weight of x; is 5, then 
the intensity for w;* is 25 pW. Of course, the value of q may be chosen 
arbitrarily. However, P(wi) leer must be greater than the saturated threshold 


of the PD. Usually, this threshold is very low, approximately 0.5 wW. Hence, 
the weight may vary over a wide range. 


For example, at xj = 1, P(x;) = 


Let us assume that there are m positive-weight input signals and n-m negative- 
weight input signals in the circuit of Figure 2(a). Then, the current after 
taking the threshold is 


a-82 
I, 2 K min{P(x,), P(w;)} - y K min{P(x,), P(w;)} * KP(T) (5) 
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Here, T = T* when +KP(T) is negative, and T = -T when it is positive. 
Considering the effect of minimum saturated threshold intensity described 
above, the threshold intensity P(T) is adjusted so that P(T*) = P(T) = 0.5 wW 
and P(T) = P(T) + 0.5 wW. This is equivalent to adding a constant maintenance 
signal H where P(H) = 0.5 wW on top of PD, to overcome the saturated threshold 
light intensity. In this case, the threshold logic output is y = 1 when 

Ig 2 0. Different logic functions can be obtained by choosing different weight 
and threshold values. 





Figure 3. Experimental Results for Input Variables Number n = 3. 
The upper three waveforms indicate x,, x2, and x,, respec- 
tively; the bottom waveform shows the output Y. 

(a) AND; (b) OR; (c) MAJORITY; (d) f(x) = x, + XX3; 
(e) f(x) = x x3 + X2(x,+xK3); (£) f(x) = xX, + KX, + X; 


Figures 3 (a) - (f) are the experimental results with n = 3. The three signals 
were supplied by three LEDs driven by three synchronized square-wave voltage 
signals produced directly by the signal generator. The peak voltage is 10 V 


and the current-limiting resistor has a value of 220 9. The weight signal y? 


and threshold T are also supplied by LEDs. By using switches and variable 
resistors, it is possible to change their signs and amplitudes. The voltage 
for the driver and the main circuit is Veg = 4.5 V. The series resistance 
values corresponding to wj = 1 and 2 are 1.2 kf and 600 0, respectively. Figure 
3(a) and (b) shows the AND and OR logic at T = 1 and 3 when wj = 1. Figure 3(c) 
shows the MAJORITY logic when wj = 1 and T = 2. When more than two of x;, x:, 
and x, are 1, the output is 1. Similarly, when w; = -1 and T = -2, we have the 
MINORITY logic. Obviously, the only way to get a more abundant logic function 
is to have different input weight values. For instance, the weight-threshold 
set (2,-1,1;1) can result in f(x) = x, + X2x;. With a set such as (1,-1,1;1), 
we have f(x) = x,x; + X2(x,+x3). With (-1,-1,1;-1) we have f(x) = x, + Xz + x3. 
These results are shown in Figures 3(d)-(f), and represent only a small 
portion of the logic functions attainable from a single threshold logic gate. 
It should be pointed out that only six out of eight of the possible 
combinations of the three input signals are shown in Figures 3(a)-(f). Two 
combinations, (0,1,0) and (1,0,1), are not included. However, this input 
signal arrangement can provide the greatest number of combinations possible. 
Therefore, Figures 3(a)-(f) can be considered to reflect the true values of 
logic operation. 
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A threshold logic gate can have multiple Ve 
binary output values such as y;, y2, -°:, [ 












































Ym. Under identical weighting, each out- —— 

put end provides a different threshold. | positive | r- Ry 

This requires the reproduction of the , fry lk Z2Ef PO; r 

weight into m parts which should not it 

interfere with each other when the a 

threshold is applied. Figure 4 shows a | —— aa Y, | Y 
logic processor with double output. The ——| #sative La | rr om Oo 
processor has two amplifier PDs with w«# —, ond —_ 
opposite output characteristics: one is T Rs 





an NPN triode T, and the other is a PNP i. 
triode T,. LED, and LED, are the output os 

ends. -V,, is the negative bias. The out- Figure 4. Schematic Diagram of 
put of LED, may be 1 only with a nega- Threshold Logic Processor With 
tive weight input and greater than Double Output Y; and Y,, Where 
i threshold value. In principle, Y, = f(x) and Y, = f(x) 


when y w,x, = T, both outputs y, and y, 


should be 1. In this case, if the weight-threshold set corresponding to y, is 
(wW,,W2,°°*,W,;T), then the weight-—threshold set corresponding to y, is 
{(—w,,-W2,°°-,-W,;-T). However, because of the cutoff property of T, and T,, also 
because it cannot make correction for a threshold simultaneously, it is only 
possible to add H on the input ends T and T° for correction. Therefore, in a 
simple manner y, and y, are opposite in nature, i.e., yy = ¥2. If the threshold 
bias maintenance light H is removed,it requires another state y, = y, = 0. In 
the following discussion, the applications of such logic devices will be 
analyzed. 


III. Discussion 


Threshold logic is not only an effective means to implement the logic 
functions in multi-variable Boolean algebra, but also can be widely used in 
digital graphic processing, mode recognition, classification of fuzzy 
functions, and neural networks. Existing optical neural network models, such 
as the Hopfield model,’ are essentially based on threshold logic. Simple 
processes such as memory connection and learning are completed by a weighted 
feedback input and threshold output between neurons.® Therefore, threshold 
logic can be directly applied in neural network research 3 a matter of fact, 
Figure 2(b) can be considered as a simple optoels nic hybrid neural 
element. It can directly add and subtract weighi« -ITATORY input and 
weighted INHIBITORY input. Hence, this neural ee: ..c can perform as a 
threshold gate. In comparison, the threshold gate in existing neural networks 
requires very complicated electronics and a computer. In addition, this neural 
element can have bipolar binary outputs of +1 and -l (see Figure 4). VLSI 
technology and free-space optical connection technology will make this type of 
neural network even more attractive. 
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Two-dimensional parallel logic processing is also a major subject in optical 
computing.® Optical threshold logic also has some potential in this area. After 
connecting an OR gate in series with the threshold logic processor for the 
double variable input shown in Figure 4, 16 Boolean logic functions can be 
implemented. In general, y, can accomplish 14 logic functions other than XOR 
[i.e., exclusive OR] and XNOR [i.e., exclusive NOR]. If threshold bias is 
removed from the circuit of Figure 4, i.e., y,; = yz = 0 is allowed to exist, 
then outputs y, and y2 represent x,-X, and X,-x,, respectively, when w, = 1, 
w, = -l, and T = 0. When an OR gate level is added in, the OR gate’s output 
y — y; + yz directly yields the XOR logic. Similarly, when w, = w, = 1 and 

T = 1, y directly provides XNOR logic. Such a logic processor can serve as a 
pixel in a two-dimensional parallel logic process. It can be directly 
controlled by an optical signal to implement real-time programmable space-— 
invariant and space-variant logic. 
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